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Kosmotropic cosolvents added to an aqueous solution promote the aggregation of hydrophobic 
solute particles, while chaotropic cosolvents act to destabilise such aggregates. We discuss the 
mechanism for these phenomena within an adapted version of the two-state Muller-Lee-Graziano 
model for water, which provides a complete description of the ternary water/cosolvent /solute system 
for small solute particles. This model contains the dominant effect of a kosmotropic substance, 
which is to enhance the formation of water structure. The consequent preferential exclusion both 
of cosolvent molecules from the solvation shell of hydrophobic particles and of these particles from 
the solution leads to a stabilisation of aggregates. By contrast, chaotropic substances disrupt the 
formation of water structure, are themselves preferentially excluded from the solution, and thereby 
contribute to solvation of hydrophobic particles. We use Monte Carlo simulations to demonstrate at 
the molecular level the preferential exclusion or binding of cosolvent molecules in the solvation shell 
of hydrophobic particles, and the consequent enhancement or suppression of aggregate formation. 
We illustrate the influence of structure-changing cosolvents on effective hydrophobic interactions by 
modelling qualitatively the kosmotropic effect of sodium chloride and the chaotropic effect of urea. 

PACS numbers: 64.75.+g, 61.20.-p, 87.10.+e, 64.70.Ja 



I. INTRODUCTION 

Most processes in living organisms are adjusted to 
function in a rather limited range of physiological condi- 
tions, and important deviations, such as high concentra- 
tions of dangerous substances, are generally expected to 
preclude life. Nevertheless, many living systems survive 
stresses of this kind and exist in hostile environments. 
One common adaptation strategy is to modify the prop- 
erties of the solvent, which is usually water, in such a way 
as to exclude the undesirable solutes from solution pj. 
Water is modified by relatively high concentrations of 
stabilising solutes (cosolvents), which remain compatible 
with the metabolism of the cell even at very high concen- 
trations (therefore they are also referred to as 'compatible 
osmolytes') These cosolvents neutralise dangerous 

solutes by decreasing their solubility and enhancing the 
formation of their aggregates. Such cosolvents are known 
as promoters of the water structure and are therefore re- 
ferred to as kosmotropes ('kosmo-trope' = order maker). 
Their stabilising function for proteins and aggregates of 
hydrophobic solute particles, and their importance for 
the osmotic balance in cells, has generated a growing in- 
terest in the physical origin of the kosmotropic effect. 
Many recent investigations have focused on the ability 
of kosmotropic cosolvents to influence solute solubility in 
aqueous media 0, E H E @ . 

Of equal interest and importance is the ability of 



chaotropic cosolvents to increase the solubility of non- 
polar solute particles in aqueous solutions 0, 0, 0, 
El EH 113 • For certain systems the solubility may be 
enhanced by several orders of mag nitude [H [H [3, 
leading to a complete destabilisation of solute aggregates 
and, in the case of protein solutions, to a complete denat- 
uration and loss of function. A significant number of sub- 
stances display this property, notably (under most condi- 
tions and for the majority of solute species) urea, which 
is frequently used in highly concentrated (c. 0.5-10 M) 
solutions as a protein denaturant. Chaotropic cosolvents 
('chao-tropic' = order-breaking) are less polar than water 
and thus break hydrogen bonds between water molecules, 
suppressing water structure formation 0, 0, |2(], l2l| . 
However, the exact physical mechanism for the changes 
in water structure, and for the consequent destabilising 
action of chaotropic cosolvents, is at present not fully 
understood. 

Because the origin of the kosmotropic and chaotropic 
effects appears to lie primarily in their influence on the 
solvent, rather than in direct interactions between co- 
solvent and solute p2l l23l |24| . a microscopic descrip- 
tion must begin with the unique properties of the aque- 
ous medium. Water molecules have the ability to form 
strong, intermolecular hydrogen bonds, and pure, liquid 
water may form extended hydrogen-bonded networks, be- 
coming highly ordered (Fig. ^| . Although the insertion 
of a hydrophobic molecule leads to a destruction of lo- 
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FIG. 1: Schematic illustration of water structure as a conse- 
quence of the formation of more or less extended networks of 
hydrogen bonds. 



cal hydrogen-bond structure, at low temperatures water 
molecules are found to rearrange in a cage-like config- 
uration around small solute particles and around high- 
curvature regions of larger ones. The orientation imposed 
on the water molecules at the interface with the particle 
results in the cage hydrogen bonds being slightly stronger 
than before, and causing a net reduction of energy dur- 
ing the insertion process [2^, |2(| |2]l 13 |2j| . However, at 
higher temperatures the free energy of the system is de- 
creased by reducing the local restructuring of water, thus 
increasing the entropy of the solvent molecules, which 
drives the aggregation of hydrophobic particles through 
a minimisation of their total surface exposed to water. 
This competition between enthalpic and entropic effects 
in the solvent is fundamental to the phenomenology of 
liquid-liquid demixing processes. The effective hydropho- 
bic interaction between small, non-polar solute particles 
is thus thought to be primarily solvent-induced, i.e. to 
be a consequence of changes in the ordering of water 
molecules rather than being controlled by direct water- 
solute interactions [30l l3ll l32l l33l l34| . In this descrip- 
tion, the primary contribution to the action of structure- 
changing cosolvents, which are generally highly soluble 
and uncharged in physiological conditions, then depends 
on their ability to alter this local ordering of liquid water. 

Kosmotropic cosolvents, such as sucrose and betaine, 
are more polar than water and act to enhance its struc- 
ture due to their ability to form hydrogen bonds @. 
For the same reason they interact with water molecules 
rather than with non-polar solute particles, which leads 
to an effective preferential exclusion from the solvation 
shell of hydrophobic molecules 0, 0, 0] , and thus to a 
stronger net repulsion between solute and solvent. In 
the presence of kosmotropic cosolvents, structural ar- 
rangement of the water-cosolvent mixture is enthalpically 
favourable compared to a cage-like organisation around 
hydrophobic solute particles. Solute molecules are thus 
pushed together to minimise their total exposed surface, 
which results in an enhancement of hydrophobic aggre- 
gation. The same process leads to a stabilisation of na- 
tive protein configurations, in spite of the fact that kos- 
motropic substances have no net charge and do not in- 
teract directly with the proteins [E 0> HI • 

Chaotropic cosolvents, such as urea in most ternary 
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FIG. 2: Schematic representation of preferential phenomena 
in a mixture of water and hydrophobic solute particles in the 
presence of chao trop ic (left) and kosmotropic (right) cosol- 
vents (after Ref. |2l1. 



systems, are less polar than water, so that their presence 
in solution leads to an energetically unfavourable disrup- 
tion of water structure. Such cosolvents are therefore ex- 
cluded from bulk water, an effect known as "preferential 
binding" to the solute particles |llj, |2J] , although it relies 
less on any direct binding of cosolvent to solute (which 
would enhance the effect) than on the fact that the cosol- 
vent molecules are pushed from the solvent into the shell 
regions of the solute. Preferential binding and exclusion 
of cosolvent molecules are depicted in Fig. [5] In the for- 
mer case, the smaller number of water molecules in con- 
tact with the surface of non-polar solute particles leads to 
a weaker effective interaction between solute and solvent, 
such that a larger mutual interface becomes favourable. 
The addition of chaotropic cosolvents to aqueous solu- 
tions therefore results in an increase in solvent-accessible 
surface area which destabilises hydrophobic aggregates, 
micelles and native protein structures 0, 0, 124. l3q | . 

The essential phenomenology of the kosmotropic and 
chaotropic effects may be explained within a solvent- 
based model founded on the concept of two physically 
distinct types of solvent state, namely ordered and dis- 
ordered water 

S II (Fig. HJ. In Sec. [iH we review 
this model, described in detail in Ref. |3q. and fo- 
cus on the adaptations which represent the structure- 
changing effects of the cosolvents. Sec. IIIII presents 
the results of Monte Carlo simulations for kosmotropic 
substances, which illustrate their stabilising effect on 
hydrophobic aggregation, the decrease in solubility of 
non-polar molecules in water-cosolvent mixtures and the 
underlying preferential exclusion of cosolvent molecules 
from the solvation shell of hydrophobic particles. We con- 
sider the inverse phenomena of the chaotropic effect in 
Sec. IIVI In Sec. we illustrate these results by compar- 
ison with available experimental data for a kosmotropic 
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and a chaotropic cosolvent, and discuss the extent of the 
symmetry between the two effects in the context of the 
model description for weakly or strongly active agents. 
Sec. IVII contains a summary and conclusions. 



II. MODEL 
A. Hydrophobic- Polar Model 

The driving force in the process of solvation and ag- 
gregation is the effective hydrophobic interaction between 
polar water and the non-polar solute |23ll30| . As noted in 
Sec. I, the origin of this interaction is the rearrangement 
of water around the solute particle. For solute surfaces of 
sufficiently high curvature to permit (partial) cage forma- 
tion, this process decreases the enthalpy due to reinforced 
hydrogen bonds between water molecules in the solvation 
shell of the hydrophobic particle in comparison to those 
in bulk water, but reduces simultaneously the number of 
degenerate states of the solvent. These physical features 
are described by the model of Muller, Lee and Graziano 
(MLG), where the energy levels and respective degenera- 
cies of water molecules are determined by the local water 
structure pslE^] . 

Because solute particles are relatively large compared 
with single water molecules, we use an adapted version of 
the MLG model in which each site contains a group of wa- 
ter molecules. The distinction between solvent molecules 
and non-polar solute particles lies in their ability to form 
hydrogen bonds. The continuous range of interaction en- 
ergies within a partially hydrogen-bonded water cluster 
may be simplified to two discrete states of predominantly 
intact or broken hydrogen bonds j3^ . The water sites in 
the coarse-grained model may then be characterised by 
two states, where an "ordered" site represents a water 
cluster with mostly intact hydrogen bonds, while a "dis- 
ordered" site contains relatively fewer intact hydrogen 
bonds (Fig. QJ. The use of a bimodal distribution in the 
adapted model is not an approximation, but a natural 
consequence of the two-state nature of the pure MLG 
model 36] . In the presence of a non-polar solute a fur- 
ther distinction is required, between "bulk" water sites 
undisturbed by the solute particles and "shell" sites in 
their vicinity whose hydrogen bonding is altered. As ex- 
plained in more detail below, this adaptation of the MLG 
model contains the essential features required to encap- 
sulate the enthalpy/entropy balance which is the basis 
for the primary phenomena of hydrophobic aggregation 
in two-component water/solute mixtures. The model has 
been used to reproduce the appearance of upper (UCST) 
and lower (LCST) critical solution temperatures and a 
closed- loop coexistence regime |40j, demonstrating that 
the origin of the hydrophobic interaction lies in the al- 
teration of water structure. 

The adapted MLG model has been extended to provide 
a minimal model for the study of chaotropic phenomena 
in ternary water/solute/cosolvent systems [36|. It was 




FIG. 3: Energy levels of a water site in the bimodal MLG 
model. The energy levels of shell water are different from bulk 
water due to breaking and rearrangement of hydrogen bonds 
in the proximity of a solute particle. The solid arrows rep- 
resent the effect of a kosmotropic cosolvent in creating more 
ordered states, and the dashed arrows the opposite effect of 
chaotropic species. 



shown that this framework yields a successful description 
of preferential binding, and of the resulting destabilisa- 
tion of aggregates of solute particles as a consequence of 
the role of chaotropic cosolvents in reducing the forma- 
tion of water structure. We begin with a brief review of 
the model to summarise its physical basis and to explain 
the modifications required for the inclusion of cosolvent 
effects. 

The microscopic origin of the energy and degener- 
acy parameters for the four different types of water 
site (ordered/disordered and shell/bulk) is discussed in 
Refs. [IE El 03- m conjunction with experimental 
and theoretical justification for the considerations in- 
volved. Here we provide only a qualitative explanation 
of their relative sizes, which ensure the competition of 
enthalpic and entropic contributions underlying the fact 
that the model describes a closed-loop aggregation re- 
gion bounded by upper and lower critical temperatures 
and densities [4£j • Considering first the energies of wa- 
ter sites (groups of molecules), the strongest hydrogen 
bonding arises not in ordered bulk water, but in the 
shell sites of hydrophobic particles due to the cage for- 
mation described in Sec. I (which can be considered as 
a result of the orientational effect of a hydrophobic sur- 
face). Both types of ordered water cluster are enthalpi- 
cally much more favourable than disordered groups of 
molecules, for which the presence of a non-polar parti- 
cle serves only to reduce hydrogen bonding still further, 
making the disordered shell energy the highest of all. As 
regards the site degeneracies, which account for the en- 
tropy of the system, cage formation permits very few 
different molecular configurations, and has a low degen- 
eracy. Water molecules in ordered clusters have signif- 
icantly reduced rotational degrees of freedom compared 
to disordered ones, with the result that the degeneracy of 
the latter is considerably higher. Finally, the degeneracy 
of a disordered shell site is still higher than that of a bulk 
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site because the relative reduction of hydrogen bonding 
due to the non-polar surface admits a higher amount of 
rotational freedom. 

Figure 13 summarises the energy levels of a water site, 
which are arranged in the sequence Eds > Edb > E a b > 
E os . Their respective degeneracies are qd s > Idb > lob > 
q os , where the states are denoted ds = disordered shell, 
db = disordered bulk, ob = ordered bulk and os = or- 
dered shell (cage conformation). We stress that these 
sequences are neither an assuption of the model nor may 
they be altered if the model is to represent a system 
exhibiting the closed-loop form of the coexistence curve 
for aggregation: only these specific sequences reproduce 
appropriately the competition of enthalpic and entropic 
contributions to the free energy. From the microscopic 
considerations of the previous paragraph these sequences 
are entirely plausible, and they have been confirmed by 
a range of experimental measurements; the exact values 
of the parameters are not important for the qualitative 
properties of the system while the order of energies and 
degeneracies is maintained. We return in Sec. V to the 
issue of departure from these sequences. 

In the calculations to follow we have used the en- 
ergy values E ds = 1.8, E d b = 1.0, E ob = -1.0 and 
E os = —2.0, which are thought to be suitably represen- 
tative for aqueous solutions, and which have been suc- 
cessful in describing a number of different types of so- 
lution 0) El EH- Their corresponding degeneracies, 
normalised to a non-degenerate ordered shell conforma- 
tion, are taken as qds = 49, qdb — 40, q b = 10 and 
q os = 1 EH ■ We remind the reader that these energy and 
degeneracy parameters are independent of temperature 
and solute density. Their relative values have been found 
to be appropriate for reproducing the primary qualitative 
features of hydrophobic interactions liol l4l] . swelling of 
biopolymers |43j, p rotein denaturation |44l |45| . and mi- 
celle formation [46[. The essential behaviour obtained 
within this model is indeed found to be rather insensitive 
to the precise parameter values, which may, however, be 
refined by comparison with experimental measurements 
to yield semi-quantitative agreement for different solu- 
tions 0) 03 Cj| • The energy scale is correlated directly 
with a temperature scale, which we define as ksT = 

We conclude this subsection by qualifying the range of 
validity of the adapted MLG model, including what is 
meant above by "small" solute particles. The model is 
not appropriate for solute species smaller than a group 
of water molecules with intercluster hydrogen-bonding, 
which necessitates a linear cluster (and solute) size of 
at least 2 water molecules. The model will also fail to 
capture the physics of large, uniform systems where the 
linear cluster size exceeds perhaps 10 water molecules, 
because in this case the majority of the water molecules 
in a "shell site" would correspond effectively to bulk 
water. Most molecular-scale treatments of water struc- 
ture in the vicinity of non-polar solute particles assume 
a spherical and atomically flat solute surface. These 
analyses 0, HE US indicate that "clathrate" models 



for the (partial) formation of hydrogen-bonded cages be- 
come inappropriate for larger curvatures, on the order 
of 1 nm, where a crossover occurs to "depiction" or 
"dangling-bond" models in which hydrogen-bond forma- 
tion is frustrated [50j ■ While the adapted MLG model be- 
comes inapplicable for such particles at the point where a 
clathrate description breaks down, we note [461 ] that most 
solute species, and in particular proteins, are atomically 
rough, in that they composed of chains and side-groups 
with high local curvatures on the order of 0.3-0.5 nm. 
As a consequence the model may remain valid for solute 
particles of considerably greater total dimension. For 
the purposes of the present analysis, we note that a sig- 
nificant quantity of the available data concerning cosol- 
vent effects has been obtained for protein solutions, and 
that the majority of this data is consistent with our ba- 
sic picture of the dominance of water structure effects in 
dictating solute solubility. However, while most proteins 
may indeed fall in the class of systems well described by a 
clathrate picture, we stress that the adapted MLG model 
cannot be expected to provide a complete description for 
large solute molecules containing many amino-acids of 
different local hydrophobicity and chemical interactions. 



B. Cosolvent Addition 

In comparison with the hydrophobic solute particles, 
cosolvent molecules are generally small and polar, and 
are therefore included directly in water sites by changing 
the number of states of these sites. A site containing wa- 
ter molecules and a cosolvent particle is referred to as a 
cosolvent site. Kosmotropic cosolvents, being more polar 
than water, increase the number of intact hydrogen bonds 
at a site . In the bimodal MLG framework their addi- 
tion to weakly hydrogen-bonded, disordered clusters may 
be considered to create ordered clusters with additional 
intact hydrogen bonds. The creation of ordered states 
from disordered ones in the presence of a kosmotropic 
cosolvent increases the degeneracy of the former at the 
expense of the latter. This feature is incorporated by 
raising the number of possible ordered states compared to 
the number of disordered states (solid arrows in Fig.[3J), 

qob,k = qob + f]bi Qdb,k = Idb — f]bi 

q os ,k = q os + r] s , qd s ,k = qds - r) s , (1) 

where k denotes the states of water clusters containing 
kosmotropic cosolvent molecules, and the total number 
of states is kept constant. 

The effect of chaotropic cosolvents on the state de- 
generacies is opposite to that of kosmotropic substances, 
and may be represented by inverting the signs in eq. Q 
Chaotropic cosolvents are less strongly polar than wa- 
ter, acting in an aqueous solution of hydrophobic parti- 
cles to reduce the extent of hydrogen bonding between 
water molecules in both shell and bulk sites [23l l5l| . 
Within the adapted MLG framework, this effect is re- 
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produced by the creation of disordered states, with addi- 
tional broken hydrogen bonds and higher enthalpy, from 
the more strongly bonded ordered clusters (dashed ar- 
rows in Fig. I3J), whence 



Qob,c — Qob — Vb, 
Qos,c Qos Vsi 



Qdb.c — Qdb + Vb, 

qds.c = Qds+Vs- (2) 



We stress three important qualitative points. First, 
the cosolvent affects only the number of intact hydrogen 
bonds, but not their strength [2j|, as a result of which, 
also in the bimodal distribution of the coarse-grained 
model, the energies of the states remain unchanged. Co- 
solvent effects are reproduced only by the changes they 
cause in the relative numbers of each type of site [eqs. (JTJ 
and J2J)]. Secondly, cosolvent effects on the relative de- 
generacy parameters are significantly stronger in the bulk 
than in the shell, because for obvious geometrical reasons 
related to the orientation of water molecules around a 
non-polar solute particle fsfil Elj , many more hydrogen- 
bonded configurations are possible in the bulk. Finally, 
as suggested in Sec. I, our general framework does not 
include the possibility of cosolvent-solute interactions, 
which have been found to be important in certain ternary 
systems, specifically those containing urea [52ll5^|. 

The illustrative calculations in Sees. Ill and IV are per- 
formed with r)b — 9.0 and r) s = 0.1 in eqs. (JIJ and J5J. In 
the analysis of Ref. [3(j these values were found to pro- 
vide a good account of the qualitative physics of urea as a 
chaotropic cosolvent; in Sec. V we will return to a more 
quantitative discussion of the role of these parameters. 
By using the same values for a hypothetical kosmotropic 
cosolvent we will demonstrate that a symmetry of model 
degeneracy parameters does not extend to a quantitative, 
or even qualitative, symmetry in all of the relevant phys- 
ical phenomena caused by structure-changing cosolvents. 
We comment briefly that the fractional value of rj s arises 
only from the normalisation convention q os = 1. 



C. Mathematical Formulation 

To formulate a description of the ternary solution 
within the framework of statistical mechanics, we begin 
by associating with each of the N sites of the system, 
labelled i, a variable n.;, which takes the values rij = 1 
if the site contains either pure water, rij = if the site 
contains a hydrophobic particle, or m = — 1 if the site 
contains a water cluster including one or more cosolvent 
molecules. This system is described by the Hamiltonian 

H[{rii}, {>i}] = ^ n »( n »+! [{E ob 8i t<Job +E db 5i^ ih )\i 



+(E 08 ~5 % ,a os + E ds S hads )(l - Ai)] 



N 



rii(ni - 1) 



[{Eob8i,(T ob „ + Edb8i,<T db a )\ 



where a = k,c denotes sites containing water and a kos- 
motropic or chaotropic cosolvent molecule. The site vari- 
able Ai is defined as the product of the nearest neigh- 
bours, \i — Tltij) an d takes the value 1 if site i is 
completely surrounded by water and cosolvent or oth- 
erwise. The first sum defines the energy of pure water 
sites and the second the energy of cosolvent sites. Be- 
cause a water site i may be in one of q different states, 
Si_ ao3 is 1 if site i is occupied by water in one of the q os 
ordered shell states and otherwise, and Si t<Tds is 1 if it is 
occupied by pure water in one of the qds disordered shell 
states and otherwise. Analogous considerations apply 
for the bulk states and for the states of cosolvent sites. 
A detailed description of the mathematical structure of 
the model is presented in Refs. an d PHI ■ 

The canonical partition function of the ternary iV-site 
system, obtained from the sum over the state configura- 
tions {ai}, is 



z n - E{„ 1 }]1 1 Z s 



i%±il (1 _ Aj) ^_i(^±l) Ai 



X Z s a 



F^(l-A 4 ) 



7 2 - 
J b,a 



\ (4) 



where Z„ 



+ 



-0E da 



for the shell (a = 



s) and bulk (a = b) states both of pure-water and of 
cosolvent sites (<r = s, a and a = b,a). 

When the number of particles may vary, a chemical 
potential is associated with the energy of particle addi- 
tion or removal. The grand canonical partition function 
of the system for variable particle number becomes 



2 = e 0^N w +0(^+A tJ ,)N aZ 



X 



(5) 



+{E OS ~5^ OB „ + E da 8 itlId . J(l - A,)], 



(3) 



N w denotes the number of pure water sites, N a the num- 
ber of cosolvent sites and N p the number of solute par- 
ticle sites, the total number being N = N w + N a + N p . 
The variable \i represents the chemical potential associ- 
ated with the addition of a water site to the system and 
A/x the chemical potential for the insertion of a cosolvent 
molecule at a water site. From the role of the two cosol- 
vent types in enhancing or disrupting water structure, at 
constant solute particle density A^i is expected to be neg- 
ative for kosmotropic agents and positive for chaotropic 



D. Monte Carlo Simulations 

We describe only briefly the methods by which the 
model may be analyzed; full technical details may be 
found in Ref. |40j . Because our primary interest is in 
the local solute and cosolvent density variations which 
demonstrate aggregation and preferential phenomena, we 
focus in particular on Monte Carlo simulations at the 
level of individual solute molecules. Classical Monte 
Carlo simulations allow the efficient calculation of ther- 
mal averages in many-particle systems with statistical 
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fluctuations > sucn as that represented by eq. (0 . We 
work on a cubic lattice of 30 x 30 x 30 sites (N = 27 000), 
with random initial particle distributions and with peri- 
odic boundary conditions to eliminate edge effects (al- 
though clearly we cannot eliminate finite-size effects). 
Each site is occupied by either a solute particle, pure 
water or a water-cosolvent mixture. The results are un- 
affected by changes in lattice size. We have implemented 
a Metropolis algorithm for sampling of the configuration 
space. After a sufficiently large number (100 000) of re- 
laxation steps, the system achieves thermal equilibrium 
and averages are taken over a further 500 000 measure- 
ments to estimate thermodynamic quantities. Coexis- 
tence lines in the p-T phase diagram are obtained from 
the transitions determined by increasing the temperature 
at fixed chemical potential p (grand canonical sampling) , 
which results in a sudden density jump at the transition 
temperature. The solute particle densities p p correspond- 
ing to these jumps yield closed-loop coexistence curves 
in the p p -T phase diagram. Despite the rather crude ap- 
proximation to a continuum system offered by the use 
of a cubic lattice, we have found quantitative agreement 
with other theoretical approaches, which is why we focus 
primarily on Monte Carlo simulations here. In Sec. V we 
will also demonstrate remarkable qualitative, and in some 
cases semi-quantitative, agreement with experiment, sug- 
gesting that the adapted MLG framework may indeed 
form a suitable basis for more sophisticated modelling. 

For the macroscopic properties of the system, such 
as aggregation and phase coexistence, the results of the 
Monte Carlo simulations may be interpreted with the aid 
of a mean-field analysis. If the densities at each site are 
approximated by their average values in the solution, p p 
for hydrophobic particles, p w for pure water and p a for 
cosolvent, where p p + p w + p a = 1, the mean value (rii) 
at site i is given by (m) = Yl,<, n i,°P° = Pw - Pa and 
analogously (nj) = J2 a n i,aP<? = Pw + Pa- The equilib- 
rium densities are obtained by minimisation of the grand 
canonical mean-field free energy per site, 

/ = (M In Z s )p w + (^ + A t i-0- 1 In Z s , a )p a 

+/3 _1 (lnZ s - In Z b )p w (p a + p w ) z 

a ~r Pw ) 

+P~ 1 (pw In p w + pa In p a + p P In p p ), (6) 

where z denotes the number of nearest neighbours of each 
site. 



III. KOSMOTROPIC EFFECT 

We begin our analysis of the kosmotropic effect by dis- 
cussing the mean-field p p -T phase diagram for different 
cosolvent concentrations, shown in Fig. As expected 
from Sec. IIA the system exhibits a closed-loop coex- 
istence curve, forming a homogeneous particle-solvent - 
cosolvent mixture below a LCST and above an UCST for 
all concentrations, whereas between these temperatures, 
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FIG. 4: Closed-loop coexistence curves for a ternary system 
of water, kosmotropic cosolvent and hydrophobic particles for 
two different cosolvent densities pk, obtained by mean-field 
calculation with degeneracy parameters r\ b = 9.0 and r\ a = 
0.1 in eq. Q. On the outside of each curve the solution 
is homogeneous, whereas on the inside it separates into two 
phases. The dashed arrow represents the heating process of a 
system with particle density po and cosolvent density pk2 = 
0.61 from temperature To in the homogeneous region to Ti in 
the two-phase region. At Ti the system is separated into two 
phases of different densities pi (nearly pure water-cosolvent 
mixture) and pi (hydrophobic aggregates). 



and for intermediate particle concentrations, a phase sep- 
aration is found into a pure solvent phase (meaning a 
water-cosolvent mixture) and an aggregated phase with 
fixed solute density. 

The hydrophobic repulsion between water and non- 
polar solute particles, developed as a result of the forma- 
tion of water structure, increases in the presence of kos- 
motropic cosolvents (this result may be shown by rewrit- 
ing the partition function (eq. |3J in terms of the par- 
ticle sites 36] ). As shown in Fig. the temperature 
and density ranges of the aggregation region indeed rise 
with increasing cosolvent concentration, illustrating the 
stabilising effect of kosmotropic agents. Further, the par- 
ticle density of the aggregate phase in the two-phase re- 
gion increases when adding cosolvent, demonstrating in 
addition the strengthening of effective hydrophobic in- 
teractions between solute particles due to the growing 
water-solute repulsion. This is illustrated by the process 
of heating a system at constant density: the dashed ar- 
row in Fig. 01 represents a solution with particle density 
Po, which is heated from temperature To in the homo- 
geneous region to a temperature between the LCST and 
the UCST. In the heterogeneous region (Ti), the solu- 
tion separates into two phases of densities p\ (almost 
pure solvent) and p 2 (hydrophobic aggregates) under the 
constraint po{V\ + V 2 ) = P\V\ + p 2 V 2 , where Vi is the 
volume occupied by phase i. An increase in cosolvent 
density leads to a higher density p 2 of hydrophobic ag- 
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FIG. 5: Closed-loop coexistence curves obtained by Monte 
Carlo simulations for an aqueous solution of hydrophobic par- 
ticles with the same parameters as in Fig. 



gregates and a lower value p± , which may be interpreted 
as a strengthening of the hydrophobic interactions gen- 
erated between solute particles. 

However, the mean-field approximation neglects fluc- 
tuations in the density, and consequently is unable to 
detect intrinsically local phenomena such as preferential 
exclusion of cosolvent particles from the solvation shell of 
a hydrophobic particle (Fig. [21 ■ Because the kosmotropic 
effect is strongly dependent upon local enhancement of 
water structure, a full analysis is required of the spatial 
density fluctuations in the system. We thus turn to the 
results of Monte Carlo simulations, which at the macro- 
scopic level show the same enhancement of hydrophobic 
aggregation (Fig. [3J as indicated in the mean-field anal- 
ysis. Below the LCST and above the UCST, one homo- 
geneous mixture is found where the hydrophobic parti- 
cles are dissolved, while between these temperatures they 
form aggregates of a given density. The LCST and the 
critical aggregate densities determined by the mean-field 
calculation do not differ strongly from the numerical re- 
sults, but the UCST lies at a temperature higher than 
that determined by Monte Carlo simulations. This is to 
be expected because mean-field calculations neglect fluc- 
tuation effects, generally overestimating both transition 
temperatures. The agreement of the mean-field result for 
small pk with the simulated value for the LCST suggests 
a predominance of local effects at low temperatures, and 
that density fluctuations between sites are small. We 
note, however, that a cosolvent concentration p k = 0.61 
is required in the mean-field approximation to produce 
an expansion of the aggregated phase similar to that ob- 
served in the Monte Carlo results for pf. = 0.25, demon- 
strating the growing importance of fluctuation effects at 
higher cosolvent concentrations. 

The p-T phase diagram of the ternary system obtained 
by Monte Carlo simulations is presented in Fig. [SJ Lines 



FIG. 6: p-T phase diagram of hydrophobic solute in an aque- 
ous cosolvent mixture for two different kosmotropic cosol- 
vent concentrations pk, obtained by Monte Carlo simulations. 
The endpoints of the finite transition lines correspond to the 
UCST and LCST. 



of finite length terminated by the UCST and LCST de- 
marcate the aggregation phase transition. An increase in 
cosolvent density leads to a larger separation of UCST 
and LCST as a direct result of the expansion of the ag- 
gregation regime (Fig. [SJ. The transition line is shifted 
to higher values of p, indicating an increased resistance 
of the system to addition of water at constant volume 
in the presence of cosolvent. This is a consequence of 
the fact that the resistance to addition of hydrophobic 
particles depends on both p and Ap: with increasing co- 
solvent concentration, Ap decreases (i.e. becomes more 
negative) and hence p must increase for the critical par- 
ticle density to remain constant. 

By inspection of the energy levels, and of the alteration 
in their relative degeneracies caused by a kosmotropic 
cosolvent (Fig.[3J), it is energetically favourable to max- 
imise the cosolvent concentration in bulk water with re- 
spect to shell water. Thus one expects a preferential 
exclusion, also known as preferential hydration, as repre- 
sented schematically in the right panel of Fig. [21 Except 
at the lowest temperatures, an increased number of water 
molecules in the solvation shell of non-polar solute par- 
ticles increases the repulsion between solute and solvent, 
resulting in a reduction of solubility, in an enhanced ef- 
fective attraction between solute particles due to the de- 
creased interface with water as they approach each other, 
and thus in aggregate stabilisation. 

Monte Carlo simulations confirm the expectation that 
the cosolvent concentration is lower in shell water than 
in bulk water. Figure shows the cosolvent shell con- 
centration compared to the overall cosolvent concentra- 
tion in the solution, where preferential exclusion of the 
kosmotropic cosolvent from the solvation shell of the hy- 
drophobic solute particles is observed. At constant chem- 
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FIG. 7: Relative concentration of kosmotropic cosolvent in the 
solvation shell of hydrophobic particles, obtained from Monte 
Carlo simulations at chemical potential /i = 2.5, exhibiting 
preferential exclusion. The effect is most pronounced for low 
cosolvent densities. 
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FIG. 8: Closed-loop miscibility curves of an aqueous solution 
of hydrophobic particles for different chaotropic cosolvent con- 
centrations p c , obtained by Monte Carlo simulations using the 
degeneracy parameters rj b = 9.0 and -q s =0.1 in eq. On 
the outside of the curves is the homogeneous, dissolved state 
while inside them is the aggregated state. As the cosolvent 
concentration grows, particle solubility increases, leading to a 
reduced coexistence regime. 



ical potential /i, a sharp drop occurs at the phase tran- 
sition temperature (Fig. [SJ , which is due to the sudden 
change in solute particle density. At very low particle 
density, a clear exclusion of the cosolvent from the solva- 
tion shell appears. However, at high particle density (for 
temperatures below the aggregation phase transition) the 
effect is only marginal, because here the number of parti- 
cle sites becomes substantial, most solvent and cosolvent 
sites are shell sites and thus the total cosolvent density is 
almost identical to the shell cosolvent density. The com- 
paratively strong fluctuations in relative cosolvent shell 
concentration at temperatures above the phase transition 
may be attributed to the very low particle density in the 
system. Small, thermal fluctuations in the number of 
shell sites occupied by cosolvent molecules result in large 
fluctuations in the cosolvent shell concentration. 

Considering next the effect of the cosolvent concentra- 
tion, the system shows a substantially stronger preferen- 
tial exclusion for low than for high cosolvent densities. 
This is again largely a consequence of considering the 
relative, as opposed to absolute, decrease in cosolvent 
concentration in the shell. Preferential exclusion is less 
pronounced for high cosolvent densities because a smaller 
proportion of cosolvent sites can contribute to the effect. 
At high temperatures, entropy, and therefore mixing, ef- 
fects become dominant and the preferential exclusion of 
cosolvent particles shows a slight decrease. 

We conclude this section by commenting also on the 
phenomenon shown in Figs. 01 and [S] that for tempera- 
tures where no aggregation is possible in the pure water- 
solute solution (below the LCST), it can be induced only 
by the addition of a cosolvent. Aggregation is driven by 
exclusion of solute from the solution in this intermediate 
regime, and is thus an indirect result of the transfor- 



mation of disordered bulk water states into energetically 
favourable ordered bulk states by the kosmotropic cosol- 
vent. Preferential hydration is a further manifestation of 
this exclusion. 



IV. CHAOTROPIC EFFECT 

A detailed study of chaotropic cosolvent effects within 
the adapted MLG model for a hydrophobic-polar mixture 
was presented in Ref. |36|. Here we provide only a brief 
review of the results of this analysis, with emphasis on 
the differences and similarities in comparison with the 
kosmotropic effect, which we have discussed at greater 
length in Sec. IIIII Figure 0D shows the coexistence curves 
demarcating the aggregated state for a range of cosolvent 
concentrations, computed by Monte Carlo simulations. 
The general features of the phase diagram are those of 
Sec. IIIII (Figs. ElandloT) . with the obvious exception of the 
shrinking of the coexistence regime as cosolvent concen- 
tration is increased. However, it appears that for the q 
and r] parameters used in the analysis (Sec. HTfe 1 ) the ef- 
fect of high chaotropic cosolvent concentrations is rather 
stronger than that obtained with the kosmotropic cosol- 
vent, an issue to which we return in Sec.0 Figure|Hlsug- 
gests that sufficiently high concentrations of chaotropic 
cosolvent may cause complete aggregate destabilisation 
at any temperature and solute particle density. 

Figure shows the corresponding fi-T phase diagram. 
As in Fig. |B| the lines representing the discontinuous 
jump in density from the low- to the high-p p phase move 
to higher /i with increasing cosolvent concentration, again 
reflecting the increased cost of adding water to the sys- 
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FIG. 9: fp-T phase diagram of hydrophobic solute in an aque- 
ous cosolvent mixture for different chaotropic cosolvent con- 
centrations p c , obtained by Monte Carlo simulation. The fi- 
nite transition lines terminate at an UCST and a LCST, which 
approach each other with increasing cosolvent concentration. 



FIG. 10: Preferential binding effect, illustrated by the relative 
concentration of chaotropic cosolvent in the solvation shell of 
hydrophobic particles, for different cosolvent concentrations 
at fixed /i = 3.0, obtained by pair approximation. 



tem in the presence of cosolvent. However, in this case 
the endpoints of the lines, which represent the UCST and 
LCST for each cosolvent concentration, become closer as 
p c increases until they approach a singularity (Fig. 
which we discuss in detail in Sec. 

Figure ^| shows as a function of temperature the co- 
solvent concentration in the shell of a solute particle nor- 
malised by the total p c . This demonstration of prefer- 
ential binding of chaotropic agents to non-polar solute 
particles was obtained from the model of eq. © within 
a pair approximation, which was shown in Ref. [36| to 
provide a semi-quantitative reproduction of the Monte 
Carlo results. As in Fig. [7| the sharp step at the phase 
transition is due to the sudden change in particle density. 
At low particle densities (below the transition) the effect 
is clear, whereas at high densities it is small because the 
majority of solvent and cosolvent sites have become shell 
sites. For low cosolvent concentrations, preferential bind- 
ing can lead to a strong increase in p c on the shell sites 
(80% in Fig. I10fl . but the relative enhancement becomes 
less significant with increasing total cosolvent concentra- 
tion. Finally, in counterpoint to the discussion at the 
end of Sec. IIHI we comment that for chaotropic agents, 
which in bulk water transform ordered states into disor- 
dered ones, preferential binding (cosolvent exclusion) is 
favoured and solute exclusion decreases. Thus a regime 
exists in the p p -T phase diagram where at fixed temper- 
ature and particle density, aggregates of solute particles 
may be made to dissolve simply by the addition of cosol- 
vent. 



V. STRONG AND WEAK COSOLVENTS 

As indicated in Sec. [H] and illustrated by the results 
of Sees. IIIII and llVl there is a qualitative symmetry be- 



tween the opposing effects of kosmotropic and chaotropic 
substances on hydrophobic aggregation of a given solute 
species. Indeed, the competition of the two effects has 
been studied for solutions of the protein HLL with the 
(denaturing) chaotrope guanidine hydrochloride and the 
(native-structure stabilising) kosmotrope betaine y| . We 
remark, however, that the same cosolvent may have a 
different effect for solutions of different types and sizes 
of solute (Sec. II), and for this reason a global charac- 
terisation of the effectiveness of any one agent may be 
misleading. 

At low cosolvent concentration, the stabilising or 
destabilising effects of the two cosolvent types on ag- 
gregate formation are rather small, and appear almost 
symmetrical in the expansion or contraction of the coex- 
istence region (Figs. and 03 Figs. Eland EJ). For weakly 
kosmotropic and chaotropic agents, by which is meant 
those causing only small enhancement or suppression of 
aggregation even at high concentrations, such qualitative 
symmetry might be expected as a general feature of the 
coexistence curves. With the same caveats concerning 
low concentrations or weak cosolvent activity, this state- 
ment is also true for preferential exclusion and binding of 
cosolvent to the solute particles, as shown by comparing 
Figs.l71andll(JI However, differing phenomena emerge for 
strongly kosmotropic and chaotropic agents. 



A. Urea as a Strong Chaotrope 

Chaotropic cosolvents act to suppress the formation of 
water structure and, because the latter is already dis- 
rupted by the presence of a non-polar solute particle, are 
more efficient for bulk than for shell sites. Within the 
adapted MLG model, a maximally chaotropic agent is 
one which causes sufficient disruption that the number 
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FIG. 11: (a) Coexistence curves for ternary water/cosol- 
vent/solute systems with different cosolvent concentrations 
p c , obtained by mean- field calculation. At a critical cosolvent 
concentration of p* = 0.81, the LCST and UCST coincide, at 
a critical temperature T* = 1.35 and critical particle density 
Pp — 0.14, where the aggregation regime is reduced to a single 
point. For still higher cosolvent concentrations the solution 
is homogeneous over the full range of temperature and par- 
ticle density, (b) Mean-field p-T phase diagram for different 
cosolvent concentrations. As p c increases, the separation of 
LCST and UCST decreases until for p c — p* c they meet at the 
temperature T£ = Tfi — T* , where the value pl ~ pu = 7.05 
determines p* v = 0.14. 



of ordered bulk states is reduced to its hypothetical min- 
imum value, which is equal to the number of ordered shell 
states; by the definition based on hydrogen-bond forma- 
tion, it is not possible to have fewer ordered bulk than or- 
dered shell states, so q ob > q os - When q ob and q os become 
similar, the entropic gain of creating ordered bulk states 
by aggregation (removal of shell sites) is largely excluded, 
and thus a strongly chaotropic agent raises the solubility 
of hydrophobic particles very significantly. With a suffi- 
cient concentration of such cosolvents, aggregation may 
be completely prevented at all temperatures and den- 
sities, and the coexistence regime vanishes at a critical 
value of cosolvent concentration. 



Such singular behaviour is well known in many aque- 
ous solutions containing urea, and indeed frequent use 
is made of the destabilising properties of this cosolvent 
(Sec. I). We stress at this point that the chaotropic ac- 
tion of urea is by no means universal: in a number of 
systems, such as the ultra-small solute methane poTIH^ . 
and proteins with variable hydrophobicity |4Sj | or signifi- 
cant interactions between chain segments and the cosol- 
vent molecules jEJ E3| , its effect may in fact be inverted. 
However, because of its ubiquity as a solubilising agent 
for most binary systems, we continue to consider urea as 
a generic example of a strong chaotrope. 

The aggregation singularity was reproduced in the 
qualitative analysis of urea as a chaotropic agent per- 
formed in Ref. 36] , and is illustrated in Fig. ^2 Because 
of the discrete nature of the Monte Carlo simulations, the 
exact parameters for the critical point cannot be found 
using this technique; these were instead obtained only 
within the mean-field approximation, and as a result 
cannot be considered to be quantitatively accurate for 
the parameters of the system under consideration. Fig- 
ure II If a) shows the shrinking of the coexistence curves 
with increasing cosolvent concentration until they vanish 
at a critical value p* = 0.81 (cf. p* « 0.55 expected from 
inspection of the Monte Carlo results in Fig. [SJl. The 
corresponding evolution of the transition line in the fi-T 
phase diagram is shown in Fig. Illf bh 

With the canonical choice of degeneracy parameters 
(Sec. EE) q ds = 49, q db = 40, q ob = 10 and q os = 1 
the values rj b = 9.0 and rj s = 0.1 required in eq. 
to reproduce this critical vanishing are indeed such that 
qob,c — 1 and q os ,c — 0.9 become very close. The level 
of "fine-tuning" of the underlying parameters required in 
the adapted MLG model reflects not a weakness of the 
framework but rather of the strongly chaotropic proper- 
ties of urea in water: the destruction of hydrogen-bonded 
networks is almost complete at high concentrations. In 
practice, concentrated urea is used to destabilise aque- 
ous solutions of many folded proteins, micelles and ag- 
gregates of hydrophobic particles. While we are unaware 
of systematic experimental studies of the effects of urea 
on the extent of the coexistence regime, aggregate desta- 
bilisation has been achieved at high concentrations in 
certain systems. At room temperature, the solubility of 
the highly hydrophobic amino-acid phenylalanin is dou- 
bled in an 8M urea solution (a solution with equal volume 
fractions of urea and water) |l6j . 



B. Sodium Chloride as a Weak Kosmotrope 

Kosmotropic cosolvents increase the extent of water 
structure formation, raising the number of ordered bulk 
and shell sites. While this process may be more efficient 
for ordered shell states, their degeneracy cannot exceed 
that of ordered bulk states. A maximally kosmotropic 
agent is one which causes sufficient structural enhance- 
ment that the number of ordered bulk states rises to the 
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FIG. 12: Coexistence curve measured in the experimentally 
accessible temperature range for N,N-Diethylmethylamine in 
water in the presence and in the absence of the kosmotropic 
cosolvent sodium chloride, reproduced from Ref. [HsL The 
LCST is reduced from 51°C in the salt-free case to — 0.6°C in 
the saturated cosolvent solution. 



point where it is equal to the number of disordered bulk 
states, q b — qdb- Singular behaviour for a strongly kos- 
motropic cosolvent (the counterpart to the vanishing of 
aggregation caused by chaotropic cosolvents) is then the 
perfect aggregation of the system with densities p% = 
and pi — 1 (Fig. i.e. extension of the coexistence 
regime across the entire phase diagram. Only at the low- 
est temperatures, where entropy effects are irrelevant, 
would the mechanism of cage formation allow solution of 
the solute particles. For a sufficiently strong kosmotropic 
effect (high concentrations of a strong agent) there is no 
reason why g Q f, should not exceed qdb, a situation which 
represents the breakdown of the sequence of energies and 
degeneracies characteristic of a system with a closed-loop 
miscibility curve (Sec.[HjY). Complete aggregation is the 
hallmark of this breakdown. 

The "symmetrical" parameters r\b = 9.0 and ry s = 0.1 
used in the analysis of Sec. IIIII give the values q b = 19 
and qdb — 31 [eq. QJ]. These are clearly not suffi- 
ciently close to represent a strongly kosmotropic cosol- 
vent, which would explain both the rather modest ex- 
pansion of the coexistence regime even with relatively 
high cosolvent concentrations (pt) and the small values 
of the upper particle density p 2 < 1 (Figs. 0] 0. One 
may thus conclude that quantitative symmetry between 
kosmotropic and chaotropic effects cannot be expected 
for any temperatures or densities, except in the event of 
very carefully chosen degeneracy parameters. 

To our knowledge, rather little data is presently avail- 
able concerning the phase diagram of hydrophobic solute 
particles in cosolvent solutions. However, Ting et al. [5^ 
have measured the lower coexistence curves for N,N- 
Diethylmethylamine both in pure water and in a mixture 
of water and the cosolvent sodium chloride (Fig. I12L 



which is usually considered to be weakly kosmotropic. 
A decrease of the LCST from 51°C to -0.6°C is ob- 
served for a solution saturated with sodium chloride, cor- 
responding to extension of the aggregation regime over a 
significantly greater temperature range. Although the 
data do not extend to sufficiently high temperatures 
to discuss the full extension of the aggregation regime 
in density (Fig. 1121 inset), the qualitative similarity of 
Figs. 01 and [5] to the experimentally determined coexis- 
tence curves (Fig. El main panel) suggests that sodium 
chloride would be appropriately classified as a "weak" 
kosmotrope in its effect on N,N-Diethylmcthylaminc. 
While its activity is significantly stronger (or its sat- 
uration concentration significantly higher) than that of 
sodium sulphate, investigated by the same authors [55| . 
at maximal concentrations it appears that the aggre- 
gation singularity is not reached and the closed-loop 
form of the miscibility curve is preserved. Here we note 
that, while the kosmotropic action of sodium chloride is 
due to its ionic nature in solution, the considerations of 
the adapted MLG framework (Sec. ITTfe) remain valid to 
model this effect. A more detailed comparison of the 
model results with the data is precluded by the fact 
that the cosolvent concentration in the saturated solu- 
tion changes with both temperature and solute density. 
Measurements of different solutions are required to char- 
acterise the dependence of aggregation enhancement on 
the molecular properties of the solute and on both solute 
and cosolvent concentrations in the system. 

Finally, preferential exclusion of kosmotropic cosol- 
vents from the immediate vicinity of hydrophobic par- 
ticles in aqueous solutions is measured only for very high 
cosolvent concentrations. Concentrations of compatible 
osmolytes, such as sucrose and betaine in the cytoplasm, 
typically reach values well in excess of 0.5M. The stabil- 
ity of the protein lactate dehydrogenase against thermal 
denaturation increases by approximately 90% in a 1M 
sucrose solution at room temperature, and by 100% for 
the same concentration of hydroxyecotine Q , classifying 
these cosolvents as being strongly kosmotropic. In the 
current model it is difficult to specify the exact cosolvent 
concentration because one site contains a group of water 
molecules, but the results are nevertheless qualitatively 
consistent with observation. We have found that, except 
for specifically chosen values of the parameters q a and 
r) a , cosolvent concentrations well in excess of 10% are re- 
quired to show a clear stabilising effect on hydrophobic 
aggregates (Fig.EJ. 



VI. SUMMARY 

We have discussed the physical mechanism underlying 
the action of kosmotropic and chaotropic cosolvents on 
hydrophobic aggregates in aqueous solutions. The aggre- 
gation phenomena arise from a balance between enthalpic 
and entropic contributions to the solvent free energy, 
which is fully contained within an adapted version of the 
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bimodal model of Muller, Lee and Graziano. This model 
considers two populations of water molecules with differ- 
ing physical properties, and generates indirectly the ef- 
fective hydrophobic interaction between solute particles. 
Such a description, based on the microscopic details of 
water structure formation around a solute particle, is ex- 
pected to be valid for small solute species, characterised 
by lengthscales below 1 nm, and for those larger solute 
particles whose solubility is dominated by the atomically 
rough (highly curved) nature of their surfaces. 

In the adapted MLG framework, one may include 
the ability of kosmotropic cosolvents to enhance and of 
chaotropic cosolvents to disrupt the structure of liquid 
water simply by altering the state degeneracies [eqs. 
and J5J] associated with the two populations characteris- 
ing this structure. As a consequence only of these changes 
in the properties of the solvent, the model reproduces all 
of the physical consequences of cosolvents in solution: 
for kosmotropes the stabilisation of aggregates, expan- 
sion of the coexistence regime and preferential cosolvent 
exclusion from the hydration shell of hydrophobic solute 
particles; for chaotropes the destabilisation of aggregates, 
a contraction of the coexistence regime and preferential 
binding to solute particles. 

The microscopic origin of these preferential effects lies 
in the energetically favourable enhancement of hydrogen- 
bonded water structure in the presence of a kosmotropic 
cosolvent, and conversely the avoidance of its disruption 
in the presence of a chaotrope. In a more macroscopic in- 
terpretation, the preferential hydration of solute particles 
by kosmotropic agents can be considered to strengthen 
the solvent-induced, effective hydrophobic interaction be- 



tween solute particles, thus stabilising their aggregates; 
similarly, the preferential binding of chaotropic agents 
to solute particles reduces the hydrophobic interaction, 
causing the particles to remain in solution. The essen- 
tial cosolvent phenomena may thus be explained purely 
by the propensity of these agents to promote or suppress 
water structure formation, and the effectiveness of a co- 
solvent is contained in a transparent way in the micro- 
scopic degeneracy parameters of the solvent states. 

Our focus on water structure is intended to extract one 
of the primary factors determining aggregation and co- 
solvent activity, and is not meant to imply that other in- 
teractions are not important. The range of solute species 
and additional contributions which are beyond the scope 
of our analysis includes extremely small solutes, large, 
atomically flat solute surfaces, solutes with differing sur- 
face hydrophobicity (such as polar side-groups) and sur- 
faces undergoing specific chemical (or other) interactions 
with the cosolvent. However, from the qualitative trends, 
and even certain quantitative details, observed in a wide 
variety of aqueous systems it appears that the adapted 
MLG model is capable of capturing the essential phe- 
nomenology of cosolvent activity on aggregation and sol- 
ubility for solutes as diverse as short hydrocarbons and 
large proteins. 
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